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Abstract. The motion of a classical pendulum in a gravitational field of strength 
g is explored. The complex trajectories as well as the real ones are determined. If g 
is taken to be imaginary, the Hamiltonian that describes the pendulum becomes VT- 
symmetric. The classical motion for this PT-symmetric Hamiltonian is examined in 
detail. The complex motion of this pendulum in the presence of an external periodic 
forcing term is also studied. 
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1. Introduction 



Many papers have been written on PT-symmetric quantum-mechanical Hamiltonians. 
In nearly all cases these Hamiltonians have had rising potentials and energies that are 
all discrete and real. An example of such a class of Hamiltonians is ^1 12 E] 

H=\p^ + x^{ixy (e>0). (1) 

In a few instances quantum-mechanical Hamiltonians having real continuous spectra 
have been examined jl]. 

For the case of PT-symmetric quantum-mechanical Hamiltonians having discrete 
spectra, the associated classical- mechanical systems have also been studied [3 El IHl Hi- 
lt has been found that the classical trajectories associated with real energies are "PT- 
symmetric (they are symmetric under reflection about the imaginary axis) and exhibit 
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remarkably interesting topological properties. These trajectories can visit multiple 
sheets of a Riemann surface. However, the classical mechanics associated with quantum 
systems having continuous spectra have not yet been examined. Thus, in this paper we 
explore the classical trajectories of the periodic potential analyzed in Ref. jlj. 

Consider a classical-mechanical system consisting of a simple pendulum with a bob 
of mass m and a string of length L in a uniform gravitational field of magnitude g (see 
Fig. Q). The gravitational potential energy of the system is defined to be zero at the 
height of the pivot point of the string. The pendulum bob swings through an angle 6. 
Therefore, the horizontal and vertical cartesian coordinates X and Y are given by 

X = Lsin^, F = -Lcos^, (2) 

which gives velocities 

X = le sin e, Y = -le cos e. (3) 

We define the potential and kinetic energies V and T as 

V = -mgL cos 6, T = \m{X:^ + F^) = \mL^e'^ . (4) 
The Hamiltonian H = T for the pendulum is therefore 

H = ^mL'^e^ - mgL cose. (5) 

Without loss of generality we set m = 1 and L = 1, and for consistency with the notation 
used in earlier papers on PT-symmetric Hamiltonians we relabel 6 —>■ x to get 

H = — g cos X, (6) 

where p = x. The classical equations of motion for this Hamiltonian read 
. dH . OH . 

x = — =p, p= — — = -gsmx. (7) 
op ox 

The Hamiltonian H for this system is a constant of the motion and thus the energy E 
is a time-independent quantity. For most of this paper we take E to be real because the 
energy levels of the corresponding quantum-mechanical system are real. 

This paper is organized as follows: In Sec. El we study the complex trajectories for 
-ff in ® when E and g are real. In Sec. Elwe study the complex trajectories for H when 
E is real but g is imaginary. In Sec. ID we examine briefly the complex chaotic motion 
that results when the energy E of the pendulum is allowed to be complex and also when 
the pendulum is subject to a periodic external driving force. 



2. Pendulum in a Real Gravitational Field 



In this section we examine the solutions to Hamilton's equations ((Zj) for the case in 
which the gravitational field strength g is real. Without loss of generality we take g = 1- 
There are three cases to consider: (1) the complex closed trajectories that result when 
\E\ < 1, (2) the complex open trajectories that occur when E > 1, and (3) the complex 
closed trajectories that we obtain when E < —1. 
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Figure 1. Configuration of a simple pendulum of mass m in a uniform gravitational 
field of strength g. The length of the string is L. The pendulum swings through an 
angle 9. We define the potential energy to be at the height of the pivot. 

2.1. Closed periodic trajectories: \E\ < 1. 

When the energy E of a simple pendulum lies in the range —1<E<1, the physically 
observed motion is that of a swinging pendulum. For the example considered in this 
subsection we choose E = 0. Since a classical turning point xq is the solution to the 
equation V{xo) = E, the classical turning points for this case satisfy cosxq = 0. Thus, 

Xq = 7i/2 + mi (n G Z). (8) 

The observed classical trajectory for this pendulum is represented by x{t) on the leal-x 
axis. This trajectory oscillates on the real axis between the turning points at Xq = ±|7r. 
However, the full array of complex solutions to the classical equations of motion is 
much richer (see Fig. In this figure we see that all classical trajectories are closed 
and periodic and lie in cells of width 2tt that fill the entire complex-x plane. As a 
consequence of Cauchy's theorem, all of the closed orbits have exactly the same period; 
namely, the period of a real physical pendulum that undergoes periodic motion on the 
real axis. This is because the period T for a particle of energy E for a Hamiltonian 
of the form H = + V{x) is given in terms of a closed complex contour integral: 
T = dx/ ^y2[E — V{x)], where C is the complex closed path of the particle. Since 
the only singularity that C encloses is the square- root branch cut joining the turning 
points, the value of T is independent of C. 

Note that the orbits in Fig. |21 that lie near the turning points are approximately 
elliptical and resemble the exactly elliptic complex trajectories of the simple harmonic 
oscillator, whose Hamiltonian is H = jp"^ + |a;^. These orbits are displayed in Fig. El 
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Figure 2. Complex classical trajectories for the simple pendulmn in These 
trajectories satisfy the equations of motion in with g ~ 1 and E ~ 0. Note that 
all trajectories are closed periodic orbits and lie in cells of width 2tt that entirely fill 
the complex plane. The usual physically observed motion of a swinging pendulum is 
represented by the line segments joining pairs of turning points on the real axis. 
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Figure 3. Complex trajectories of the simple harmonic oscillator, whose Hamiltonian 
is H — + ia;^ with E = 1. These orbits are exact ellipses. These ellipses closely 
approximate the complex orbits in the vicinity of the turning points in Fig. |2] 



2.2. Open trajectories: E > 1. 

When the energy E of a. simple pendulum is sufficiently large (here, E > 1), the 
physically observed motion is that of a rotating pendulum. The classical turning points 
for this case are no longer on the real axis because there are no real solutions to 
cosxq = E with E > 1. The rotating pendulum has no real turning points because 
its real trajectory never turns back. To find the complex turning points, we replace Xq 
hj a + ib and get two equations: 

sin a sinh 6 = 0, cos a cosh b = —E. (9) 
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For the example considered in this subsection we choose E = coshl, and thus the 
turning points are situated at a = {2k + l)7r, where G Z, and 6 = ±1. 

The complex trajectories for this case are shown in Fig.|3] Note that the trajectories 
fill the entire complex plane and that there are no closed orbits. The turning points 
have a strong infiuence on the path of the complex pendulum bob. The trajectories in 
the upper-half complex-x plane, for example, go just below and veer around the turning 
points above the real axis. The turning points seem to attract the trajectories in much 
the same way as was found in Ref. [0]. 



Im(x) 




Re(x) 



Figure 4. Trajectories in the complex- a: plane for the rotating pendulum. For this 
case the pendulum has energy E = coshl and the gravitational field strength has 
the value g = I. Note that the trajectory on the real axis corresponds to a physical 
pendulum undergoing continuous circular motion. All of the complex trajectories are 
open. The trajectories that terminate at turning points run vertically off to complex 
infinity. One such trajectory that starts at xq = tt + i is shown. 



Note that the trajectories that terminate at turning points run vertically off to 
complex infinity. This feature is also a property of the Hamiltonian H = + ix^, 
which was studied in detail in Ref. [2]. The complex trajectories for this Hamiltonian 
are shown in Fig. El There are three turning points, and the trajectory beginning at the 
turning point on the imaginary axis goes up the imaginary axis to infinity. 

One may ask how long it takes for a particle beginning at a turning point, say, the 
turning point at it + i, to reach infinity. The time T for this motion is represented by 
the following integral: 

T = ^ , . (10) 

V2 A=i+7r VE + cosx 

A numerical evaluation of this integral gives T = 1.97536 ■ ■ ■. Thus, the particle reaches 

the point at complex infinity in a finite time. 




2.3. Closed trajectories: E < —1 

There are no physically observable trajectories for this choice of the energy. However, 
the complex trajectories in this case can be calculated numerically (see Fig. ^ and we 
find that they are all closed orbits that lie in periodic cells of width 2tt. The turning 
points come in complex-conjugate pairs. In Fig. IHlwe have chosen E = — cosh 1 and the 
turning points are situated at Xq = a + ib, where a = 2kTi {k G Z) and b = ±1. 



3. Pendulum in an Imaginary Gravitational Field 

If we take the gravitational field strength in Q to be imaginary, the Hamiltonian 
becomes VT symmetric, where V reflection consists of x — ^rni — x. Setting g = i, the 
Hamiltonian becomes 

H = p'^/2-icosx (11) 

and the turning points are located a.t xq = a + ib, where 

cosacosh6 = 0, — sin a sinh 6 = i?. (12) 

There are now two cases to consider, E > and E < 0. In Fig. Q we take 
E = sinh 1 and find turning points at Xq = (n + ^) vr + i(— 1)" (n G Z) and in Fig. |Hlwe 
take E = — sinh 1 and find turning points at Xq = (n + vr + i(— 1)"+-*^ (n G Z). Like 
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Figure 6. Complex trajectories for tlie pendulum in © with g = I and E — — coshl. 
The turning points come in complex-conjugate pairs and the trajectories are all closed 
orbits. 



the complex trajectories in Fig. HI the trajectories in these two figures are not closed 
and trajectories beginning at turning points run vertically off to infinity. In Fig. [7| one 
such trajectory is shown. This trajectory begins at |7r + i. The time T for a particle to 
travel from this turning point to infinity is given by 

1 r°°+3-/2 dx , , 

T = ^ , (13) 

V2 Jx=i+37r/2 V-E + 2 COSX 

After a simple change of variables this integral becomes 

T = ^ / ^ = 4K(-l/e2) = 1.84549 ■ • • , (14) 

v2 Js=i vsinh s — sinh 1 \/e 

where K is the complete elliptic function. Thus, the classical particle reaches complex 

infinity in finite time. 



4. Further Considerations 



There are a number of interesting ways to extend and generalize the investigations in 
the previous two sections. 

In Figs, m and IHl the energy for the pendulum in an imaginary gravitational field was 
taken to be real. This is because a quantum-mechanical "PT-symmetric Hamiltonian 
has the virtue of possessing real energy levels. However, one might ask what happens 
if one takes the energy to be complex. In this case the complex trajectories are not 
spatially or temporally periodic, even though the turning points are periodic and are 
separated horizontally by 27r. The behavior of the trajectories is illustrated in Fig. IHl 
where we have taken the energy to be = z. 

Another way to extend these investigations is to study what happens when there 
is an external periodic driving force of the form esin(cA;t) jHlElElE]- We have taken 
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Figure 8. Complex trajectories for the pendulum in an imaginary gravitational field 
of strength g — i and negative real energy E = — sinhl. Trajectories that begin at 
the turning points run off to complex infinity in the positive and negative imaginary 
directions as in Figs. 01 and[7| 



e = 0.2 and u = 0.1 and in Figs. ^1 ^2 ^ind ^1 we have examined the complex 
trajectories up to times t = 100, t = 650, and t = 1000. The starting energy for 
these figures is taken to be and the value of is 1. The initial value of x is taken 
to be ^TT + 0.1,which would lead to a complex trajectory if there were no driving 
term. We see that for short times the trajectory is noisy, but seems to be confined 
to the central cell in Fig. El However, after a sufficiently long time the trajectory 
wanders into adjacent cells, as can be seen in Figs. ^2 and El This result suggests it 
might be productive to investigate complex trajectories of dynamical systems whose real 
trajectories exhibit chaos, especially in the PT-symmetric case. One could begin with 
the complex trajectories associated with real homoclinic chaos, as in the present case. 
In addition, the rich behavior seen in the complex spatial trajectories of PT-symmetric 
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Figure 9. Complex trajectories for a pendulum in an imaginary gravitational field of 
strength g = i and complex energy E — i. In this case the trajectories are not closed 
and not spatially periodic. 

systems may be expected to be the projection onto the two-dimensional complex plane 
of interesting behavior in four- dimensional complex phase space. For example, the 
phase-space trajectories of PT-symmetric complex systems whose real parts are chaotic 
might exhibit a complex Hamiltonian analog of Arnold diffusion, or some other kind of 
Hamiltonian diffusion. 
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Figure 10. Complex trajectories for a pendulum in a gravitational field of strength 
g — \ and subject to an external periodic driving force. Here, the driving force has 
the form 0.2 sin(O.lt) and the initial condition is taken to be x = ^tt + 0.1. Ordinarily, 
such an initial condition would lead to a complex trajectory like that in Fig. 12 For 
short time (up to t = 100 in this figure) the trajectory remains confined to the central 
cell in Fig. H 



We thank D. C. Brody for many useful discussions. CMB thanks the Theoretical Physics Group at 



Complex Trajectories of a Simple Pendulum 

lm(x) 
7.5 



10 




5 10 15 20 25 30 35 



Re(x) 



Figure 11. Same as in Fig. 1101 but now the trajectory is allowed to continue until 
t — 650. Eventually, the trajectory excapes from its cell and begins to wander into 
nearby cells. It gains complex energy and it begins to resemble the trajectory in Fig.|51 
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Figure 12. Same as in Fig. but now the time runs up to i = 1000. Note that the 
complex trajectory undergoes extreme fluctuations. 
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